Probing Fermionic Condensates by Fast-Sweep Projection onto Feshbach Molecules 
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Fast sweep projection onto Feshbach molecules has been widely used as a probe of fermionic 
condensates. By determining the exact dynamics of a pair of atoms in time varying magnetic fields, 
we calculate the number of condensed and noncondensed molecules created after fast magnetic field 
sweeps from the BCS to the BEC side of the resonances in 40 K and 6 Li, for different sweep rates and 
a range of initial and final fields. We discuss the relation between the initial fermionic condensate 
fraction and the molecular condensate fraction measured after the sweep. 

PACS numbers: 03.75.Ss, 03.75.Kk, 05.30.Fk 



Ever since the first realisations of fermionic conden- 
sation in two-component gases of 40 K [l[ and 6 Li Q 
the question how these fragile states can be accurately 
probed has been of great experimental importance. Since 
most techniques rely on the time-of- flight expansion, dur- 
ing which the Cooper pairs dissociate, indirect meth- 
ods have been proposed. One of the most common is 
the so called fast-sweep projection P, which consists 
in a rapid magnetic field ramp from the BCS to the 
BEC side of the resonance. The key point is to make 
the field variation fast with respect to the many-body 
time scales thus preventing pair collisions from changing 
the pairs' momentum distribution and re-thcrmalisation. 
With this condition, the condensed molecules measured 
after the sweep originate only from the Cooper pairs, 
already present at the initial time, whereas the noncon- 
densed molecules are converted only from the noncon- 
densed atoms. The expectation is that this allows the 
information about the initial state to be easily unrav- 
elled from the measurements of molecules created dur- 
ing the sweep. The aim of this Letter is to describe the 
fast sweep regime of fermion dynamics and to determine 
whether and how the information about the initial con- 
densate can be extracted from the production efficiency 
and the distribution of molecules after the sweep. 

Despite its experimental relevance and fundamental in- 
terest, understanding of the dynamics of a cold Fermi gas 
following variations in the interaction strength is still rel- 
atively limited. The early analyses of molecular produc- 
tion were based on overlapping the initial pair state with 
the final molecular wave function This corresponds 
to an abrupt jump of the magnetic field strength and thus 
cannot account for the sweep rate dependence. The anal- 
ysis of time-dependent processes has so far concentrated 
on the mean-field dynamics which by definition ig- 
nores the noncondensed molecules. Thus, such methods 
are inadequate for sweeps where large numbers of non- 
condensed molecules were reported to be created during 



the ramp Q . Finally, in the dependence of the molec- 
ular production on the sweep rate was estimated analyti- 
cally without determining the actual dynamics. This was 
done by overlapping the initial state with a sweep-rate- 
dependent "final state" molecular wave function. Clearly, 
more rigorous calculations, which would determine the 
dynamics of the gas in time- varying fields and allow us 
to capture both condensed and noncondensed molecules 
created during this process, are still missing. 

Making use of a systematic cumulant expansion jjjj], 
one can obtain the dynamic equations for the density 
matrix and the pair function of the gas. On this basis, 
the number of molecules beyond the mean-field level, i.e. 
accounting for noncondensed pairs, can be calculated. In 
this Letter we numerically solve these equations for a 
limited case of fast sweeps, during which pair collisions 
are insignificant. In such a case, the dynamical aspects 
are captured by the two-body time evolution of a single 
atom pair in the spirit of the original rationale behind 
the fast sweep projection technique. We determine the 
two-body time evolution exactly by numerically solving 
the time-dependent Schrodinger equation for a pair of 
atoms. We start from the single-channel Hamiltonian 
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and a ks (hereafter k, k', and 



q denote wave vectors in three dimensions). In the 
following we use a finite-range potential of the form 
14k' = Vb(£) X (k)x(k'). Here x (k) = cxp[-(k ( 7 bg ) 2 /2] 
and the parameters Vq(B) and <7b g are chosen to recover 
the magnetic field dependence of the scattering length 
a(B) = o,b g (l — B A B Ba ) and of the highest vibrational 
bound state energy beyond the universal regime (details 
in [1, 0, [n|). The density of molecules at time i 11 1 is 
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where </>b(-B,k) is the Feshbach molecule wave function 
(determined by solving the stationary Schrodingcr equa- 
tion for a pair of atoms with the finite-range potential 
Vkk') an d the time dependence of the field is given by 
B(t) = B + B(t - t ). To calculate © with the time- 
dependent Hamiltonian JT]) we use the cumulant expan- 
sion 0]. In the presence of a condensate, higher order 
correlation functions are still large, whereas higher order 
cumulants, which include the equivalent order of inter- 
actions around the interaction-free evolution, are small - 
at least at short times. The expansion is thus controlled 
by how fast the resonance is crossed. Having expressed 
@ in terms of cumulants, we can close the hierarchy 
of their dynamical equations by neglecting all fifth- and 
higher-order terms (the second-order expansion). Note, 
however, that since here (aj^) = (os ks ) = 0, we need to 
keep only the second-order cumulants (which here are 
equivalent to the correlation functions: the density ma- 
trix f(k,£) = \ X) s ( a iLW a ksW) an d the pair function 
$(k, t) = (a_kj.(*) a kT(*))) an d the fourth-order ones of 
the form (Fii^i^-Fk) (where Fi indicate any fermionic 
creation and annihilation operators). Next, by moving 
to the interaction picture, it is possible 0] to formally 
solve for all fourth-order cumulants and to obtain close 
equations for r(k, t) and $(k, t). This results in addi- 
tional non-Markovian collision terms containing products 
of r(k, t) at different momenta and times in addition 
to terms already present in the mean-field equations [Bj]. 
Consequently, the molecular density becomes 



density of condensed molecules becomes 
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The two-body transition amplitude is defined as 
T(k,M ) - (4> h (t)\U 2 B(t,t )\k), where U 2B {t,t ) is a 
two-body evolution operator. The first term in ^ de- 
scribes condensed molecules and the remaining terms the 
noncondensed molecules created after time t. Second- 
order cumulant expansion includes two-particle collisions 
but neglects higher order collision terms. Certainly, for 
very slow ramps, for which the re-thcrmalisation and 
multiple collisions are important, a higher order expan- 
sion may be needed. However, already the solution of 
full non-Markovian second order cumulant equations is 
numerically challenging and is beyond the scope of this 
Letter (although a suitable parallelisation of the compu- 
tation should allow slower field variations to be studied). 
Here, our intention is to focus on the fast sweep limit. 
On such short time scales the evolution of the pair func- 
tion <f> is principally captured by the two-body evolution 
operator U 2 b, i- e &(t) = U 2 B(t,to)Q(to), an d the third 
(collision) term in equation ((3]) is small. In this limit the 



n c m (t) 
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(4) 



where $(to) is the initial pair function. Density of non- 
condensed molecules is 

< c = 2 fdk' |T(k', t, t )\ 2 [ dkf (k, t )f (k-2k', t ), (5) 



where T(io) is the initial density matrix. We 
determine U2B{t,to) exactly by numerically solving 
ifr-§tU 2 B(t,to) = H 2B (t)U 2B (t,t ), where H 2B is equal 
to (p} for a pair of atoms. We are then in a position to 
evaluate n c m and n^ c from j4| and j5]), respectively. Since 
our main focus is on establishing the relation between the 
final state after a fast sweep and the initial state, we take 
the simple mean- field thermodynamic initial conditions. 

We first analyse the dependence of n c m , n^f and n m = 
n m+ n m 0n the final field (Fig. [J). For finite-rate sweeps 
(unlike for the abrupt jump) the number of produced 
molecules depends on the final field only if it is in the 
region close to the resonance - it saturates further away 
from it; and the faster the sweep the further from the 
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FIG. 1: (color online). Right: The scattering length (dashed- 
doted lines) and the bound state energy, E^(B), of the 40 K 
(top) and the 6 Li (bottom) Feshbach molecule 4>b(B) in the 
vicinity of the 202 G ( 40 K) and the 834 G ( 6 Li) resonances 
from two- (solid lines) and single-channel approaches (dashed 
lines). Left: Total (top curves, blue), condensed (middle 
curves, red) and noncondensed (bottom curves, black) molec- 
ular densities versus the final field, Ba n , for the initial field 
B ini - B = 0.12 G for 40 K and B ini - B = 4.62 G for 6 Li 
and slew rates (1/B): 10 jus/G ( 40 K), 0.01 /xs/G ( 6 Li) (solid 
lines); 1 fJs/G ( 40 K), 0.003 /is/G ( 6 Li) (dashed lines); and 
abrupt jump (dotted lines). The initial atomic densities are 
1.5 ( 40 K) and 2.9 ( 6 Li) 10 13 cm -3 . For 40 K the parameters 
are from [l[ and for 6 Li are motivated by 0, 03] ■ 
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resonance this occurs. This supports the picture used in 
Q that sufficiently far from the resonance Eb is so large 
and interactions so small that the molecular state adia- 
batically follows the ramp (no molecules are created or 
dissociated during this part of the dynamics) For fi- 
nal fields near the resonance the molecular production is 
independent of the sweep rate and close to the one after 
jump. There is, however, an intermediate region of fields 
for which the molecular production is not yet saturated, 
but it depends on the rate. The presence of this region 
and also the exact position of the crossover to adiabatic 
dynamics is expected to be the main source of quanti- 
tative discrepancy between the dynamics presented here 
and analytical estimates from Q for condensed pairs (ad- 
ditional differences arise for noncondensed ones - these 
are discussed later). Note, that there is a qualitative 
difference between the 202 G resonance in 40 K and the 
834 G resonance in 6 Li [l^, El (see Fig [J). For 6 Li and 
slew rate l/B = 0.003/zs/G used in molecular 
production saturates around 250 G below Bq, where the 
single and two-channel predictions already differ (there 
is around 15% difference in Ef> between the two). This 
suggests that if we were to study sweeps even faster than 
those in [121 ] an extension to two-channels would be nec- 
essary. In contrast, for 40 K Eb from the single- and 
two-channel approaches coincide in the relevant region 
B — Bq > — 10 G and are close for any field. Thus, the 
single channel theory is sufficient for analysis of molecu- 
lar production in 40 K for any sweep rate. 
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FIG. 2: (color online). Left: Molecular condensate, 2n^/n, 
(solid lines), non-condensate, 2n^ /n, (dashed lines) and to- 
tal, 2n m /n, (doted lines) conversion efficiencies versus the 
initial field, Bi n i, for the final field Bfi n — Bo — — 10 G 
and two slew rates; and for an abrupt jump, which 

produces virtually no noncondensed molecules (bottom solid 
curve, green). Top solid curve (black) shows the initial (at 
Si n i) condensate fraction, 2n c /n. Right: Final molecular con- 
densate fraction, n^/rim, for the same parameters versus Bmi 
(top), and versus the initial 2n c /n (bottom). Bottom solid 
curves (black) show the initial 2n c /n . Calculations are for 
40 K 202 G resonance and atomic density of 1.5 10 13 cm -3 . 



The dependence of n c m on the initial field (Fig. [2]) 
is related to the number of Cooper pairs, n c (B) = 
j dk|$(k, B)\ 2 , at that field, but it is different for differ- 
ent sweep rates. As expected, the number of condensed 
molecules, n^, created during a linear sweep is always 
smaller than the number of initial Cooper pairs but larger 
than what one would obtain after an abrupt jump. The 
number of noncondensed molecules, n^f , increases as the 
initial field is shifted deeper to the BCS side. This can be 
understood by noting that, in the fast sweep regime, the 
noncondensed pairs are created only from atoms which 
were initially out of the condensate and the number of 
these atoms increases towards the BCS side. 

Finally, the dependence on the sweep rate is shown 
in Fig. [3] and compared with [|| and, for n m , with the 
Landau-Zener (LZ) theory from Since Cooper pairs 
convert to molecules more efficiently than free atoms (as 
pointed out in Q), LZ, which does not distinguish be- 
tween the two, is unlikely to give an accurate account of 
molecular production from a partially condensed Fermi 
gas. Indeed, LZ results differ substantially from the re- 
sults of our calculations (see Fig. [3]). The differences with 
Q are pronounced especially for the production efficiency 
of noncondensed molecules, which in [f| has a linear de- 
pendence on the slew rate with a coefficient dependent on 
the density but independent of the initial field. Here, it 
is instructive to compare our dynamics with an approach 
based on the asymptotic dissociation spectrum given by 
(B12) in 15|. It can be shown [l5| that in the asymptotic 
-co and t — > +oo the transition probability 
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this form instead of the exact one in ([5]) gives an approx- 
imate number of noncondensed molecules, good for B ln i 
deeply in the BCS regime. An interesting case is that of 
Bini = Bq for which obtained using the exact tran- 
sition probability, |T(k, t, to)\ , turns out to be a half of 
nj^ c obtained using the asymptotic spectrum. On the ba- 
sis of numerical evidence we can thus provide a simple 
formula for the number of noncondensed molecules for 
B m i = Bq as a function of the slew rate 



2TTh 2 m\B\ 



dk 1 cxp 



4obgAB k , 3 
3 h 2 mB 



y"dkf(M )f(k-2k',; ) (6) 



(x symbols in Fig. [3|), whose agreement with the exact 
one is very good indeed. The dependence on 1/|J5| is a 
product of a linear term, dominant for very fast sweeps, 
and an exponential correction. The coefficient of this de- 
pendence is determined by the overlap of the transition 
probability and the initial density matrix at different mo- 
menta, and thus varies with the initial field (see Fig. [2]). 

As shown in Fig. [2] and [3] (Inset) , the molecular con- 
densate fraction, n c m jn m , is always much larger than the 
initial one, 2n c /n, but also it decreases to some extent 
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FIG. 3: (color online). Conversion efficiency for condensed, 
2n^,/n, noncondensed, 2n^ c /n, and total, 2n m /n, molecules 
versus the dimensionless inverse sweep rate X for the 

40 K 202 G resonance and atomic density of 1.5 10 13 cm~ 3 for 
Bini = Bo and Ba n — Bo = -10 G. For comparison, formu- 
las from [(| (dotted lines: total, condensed and noncondensed 
from top to bottom); and for 2n m /n the Landau-Zener pre- 
diction from 0] (dashed-dotted line) are shown, x symbols 
indicate the asymptotic approximation as discussed in the 
text. Inset: Final molecular condensate fraction, n™/n m , 
(solid line) for parameters as in the main Fig. For compari- 
son, n^/rim from [|| (dotted line) and the initial condensate 
fraction, 2n c /n, (flat dashed line) are plotted. 

with the slew rate. However, for slower sweeps the many- 
body effects become more important, which would again 
lead to an increase of the condensate fraction by for- 
mation of new condensed molecules from initially non- 
condensed atoms. This effect, together with three- and 
four-body collision losses which are density dependent 
and thus larger for condensed pairs occupying a smaller 
region in the centre of the trap [l6j|, is a possible expla- 
nation for the sweep rate independent measurement of 
the condensate fraction in 40 K (Fig. 5 in [l[). Note, that 
only the fastest sweeps in this Figure are expected to be 
in the fully two-body regime. Due to the large interchan- 
nel coupling it has proven to be much harder to realise the 
fast sweep limit for 6 Li. For the fastest sweeps reported 
[l2| of 0.003 fis/G (dashed lines in the left bottom panel 
of Fig. [I}, atomic density n=2.9 10 13 cm -3 and initial 
field around 4G above the resonance, our method pre- 
dicts around 0.49 for 2n^/n and about 0.03 for 2n^/n. 
To conclude, wc have shown that even in the limit of 



fast sweeps, where the dynamics of the gas is predomi- 
nantly governed by the two-body evolution, the number 
of condensed molecules as well as the molecular conden- 
sate fraction after the sweep depend on both the sweep 
rate and the initial state. This makes it difficult to un- 
ravel the initial condensate fraction (and the initial num- 
ber of Cooper pairs) from the measurement of the final 
molecular condensate without detailed dynamical calcu- 
lations. However, although in the fast sweep limit the 
measured molecular condensate fraction always overes- 
timates the initial condensate fraction, the presence of 
the molecular condensate implies the existence of the 
fermionic condensate before the sweep. 
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